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Abstract 

A statistical mechanics model for a faceted crystal is the 3D Ising 
model at zero temperature. It is assumed that in one octant all sites 
are occupied by atoms, the remaining ones being empty. Allowed atom 
configurations are such that they can be obtained from the filled octant 
through successive removals of atoms with breaking of precisely three 
bonds. If V denotes the number of atoms removed, then the grand 
canonical Boltzmann weight is g^, < g < 1. As shown by Cerf and 
Kenyon, in the limit q ^ 1 a. deterministic shape is attained, which has 
the three facets (100), (010), (001), and a rounded piece interpolating 
between them. We analyse the step statistics as g — > 1. In the rounded 
piece it is given by a determinantal process based on the discrete 
sine-kernel. Exactly at the facet edge, the steps have more space to 
meander. Their statistics is again determinantal, but this time based 
on the Airy-kernel. In particular, the border step is well approximated 
by the Airy process, which has been obtained previously in the context 
of growth models. Our results are based on the asymptotic analysis 
for space-time inhomogeneous transfer matrices. 

1 Introduction 

As a very common phenomenon, crystals are faceted at sufficiently low tem- 
peratures with facets joined through rounded pieces. Of course, on the atomic 
scale the crystal surface must be stepped. These steps meander through ther- 
mal fluctuations. On a facet the steps are regularly arranged except for small 
errors, whereas on a rounded piece the steps have more freedom to fluctu- 
ate. Our aim is to understand the precise step statistics, where the step 
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Figure 1: Crystal corner viewed from the {111) -direction for q = 0.98. 



bordering the crystal facet is of particular interest. To gain some insight 
we will study a simplified statistical mechanics model of a cubic crystal. Its 
equilibrium shape has three facets, each consisting of a part of one of the 
coordinate planes. The facets do not touch each other and there is an inter- 
polating rounded piece, see Figure [U For this model the step statistics will 
be analyzed in great detail. 

Let us first explain our model for the corner of a crystal. The crystal is 
assumed to be simple cubic with lattice Z^. We use lattice gas language and 
associate to each site x G Z"^, Z+ = {0,1,2,...}, the occupation variable 
Ux = 0,1 with 1 standing for site x occupied by an atom and for site x 
empty. Up to a chemical potential the binding energy of the configuration n 
is 

H{n) = J J2 ^>0. (1.1) 

\x-y\=l 

We consider very low temperatures, meaning that all allowed configurations 
have the same energy, i.e. the same number of broken bonds. To define 
properly, we introduce the reference configuration n^^^ in which only the 
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octant Zi is occupied, 



1 for X G 



1 for X G \ Z^ . 

n is an allowed configuration if for a sufficiently large box A one has 

= nf for all x G Z^ \ A and H{n) - Hin'"^) = 0. (1.3) 

The set of allowed configurations is denoted by Q. By construction Q is 
countable. To favor a crystal corner, we introduce the fugacity g, < g < 1, 
and assign to each n E fl the weight 

(1.4) 

where V{n) is the number of atoms removed from rf^^, i.e. 

Vin)= J](l-n.). (1.5) 



A configuration n G can uniquely be represented by a height function 
h over Z^. For the column at G Z^, all sites below h{i,j), excluding 
h{i,j), are empty and all sites above h{i,j) are filled, n E fl if and only if 

h{i + lj) <h{ij), h{ij + 1) < h{i,j), h{i, j) ^ ioi {i, j) ^ oo. 

(1.6) 

By abuse of notation, the set of height functions satisfying ()1.6|) is also de- 
noted by Q. For h E Q let V{h) be the volume in Zi^ below h. Then the 
weight for the height h is q^^^\ 

There is an alternative way to describe configurations n E Q, which we 
just mention for completeness, but will not use later on. One builds the 
crystal out of unit cubes and projects its surface along the (lll)-direction, 
which results in a tiling of the plane with lozenges (rhombi) oriented 
along 0, 27r/3, and 4n/3. With the orientation of Figure El there are three 
sectors of the plane corresponding to the polar angle 9 with —n/G < 9 < 7r/2, 
7^/2 < 9 < 77r/6, 77r/6 < 9 < ll7r/6. n G if and only if the tiling in each 
sector becomes regular sufficiently far away from the origin. 

Instead of tilings, if preferred, one can also think of covering the dual 
hexagonal lattice by dimers such that every site is covered. In computer 
science this is called perfect matching. Equivalently, to have a more statistical 
mechanics fiavor, one can consider the fully frustrated antiferromagnetic Ising 
model on a triangular lattice, i.e for an allowed spin configuration every 
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Figure 2: The (111) projection of a configuration n In each of the three 
sectors the tiling becomes regular far away from the origin. 



triangle must have exactly two spins of the same sign. Erasing all bonds 
connecting equal sign spins yields a lozenge tiling, and viceversa. 

Following the conventional pattern we should now state our main results. 
This would mean to bring in a lot of additional notation making our intro- 
duction unwieldy. Thus we sketch only the general goals, explain our results, 
and connect to previous studies of the model. The precise theorems will be 
given in the respective sections. 

The step statistics is studied in the limit g — > 1. Thus it is convenient to 

set 

g=l-i, T^oo. (1.7) 
Let Ht denote the random height function distributed according to 

^exp[ln(l-l)F(/i)] (1.8) 

relative to the counting measure on Q, the normalizing partition function. 
For large T the heights are O (T). Thus one expects a limit shape on scale 
T. In fact, as proved in [211121, 

lim ^hriluT], [vT]) = h^^iu, v) (1.9) 
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Figure 3: Zoom to the facet edge in Figure^ 



in probability. Here {u,v) G and [ ] denotes the integer part. Let 
V = {{u,v) e ]R^,e~"/^ + e""/^ > 1}. On V, h^a. is strictly decreasing in 
both coordinates and /ima > 0, whereas /ima = on \ V. The analytic 
form of /ima is given in Section IHl If r denotes the distance to &D = {{u, v) G 
K+,e~"/^ + e"""/^ = 1}, it follows that /ima vanishes as r^/^. This is the 
Pokrovsky-Talapov law. 

Our interest here is to zoom to the atomic scale. One possibility is 
to consider a macroscopic point {u, v) ^ T> and the local height statistics 
{hT{[uT] + i, [vT] e 1?}. In the limit T ^ oo, locally the height 

profile is planar and one expects that the height statistics corresponds to a 
random tiling of the plane with the three types of lozenges from Figure |2l 
such that the relative fraction of lozenges yields the average slope V/imal^, v). 
This property will be proven in Section H] and we refer to it as local equilib- 
rium: as T — s> oo, locally one sees an infinite volume translation invariant, 
spatially ergodic Gibbs measure for the lozenges with their chemical poten- 
tials determined through V/ima(w,f). 

An even more intriguing issue is to zoom to the facet edge, which means 
to take (m, v) G dV, see Figure El Since the step density vanishes at dfD, 
typically there will be only a few steps in focus. Thus it is more natural to 
consider directly the crystal step bordering the facet. By symmetry we can 
choose the border step lying in the 2 — 3 plane. Then the border step is given 
as the graph of the function 



t^bT{t) = hT{0,t), teZ 



(1.10) 
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From ()1.6|) we have + 1) < ^t(^) and linif^oo &t(^) = 0. For large T, bj- 
is C (T), and there is a hmiting shape according to 

hm T-'bT{[TT]) = bUr), (1.11) 

T — >oo 

where 

6^(r) = -21n(l r > 0. (1.12) 

()1.12|) tells us only the rough location of the step. For the step statistics 
the relevant quantity is the size of the fluctuations of 6r([rT]) — Tb^{T). 
As will be shown they are of order T^/^ which is very different from steps 
inside the rounded piece of the crystal which are allowed to fluctuate only 
as InT. On a more refined level one would like to understand correlations, 
e.g. the joint height statistics at two points t and t' . They have a systematic 
part corresponding to boo{T). Relative to it the correlation length along the 
border step scales as T^/'^, which reflects that on short distances the border 
step looks like a Brownian motion. Thus boo has to be expanded including 
the curvature term and the correct scaling for the border step is 

T-'/'{bTi[rT + sT'/'])-{bUr)T + b'^iT)sT'/' + lb'^^^^^^ 

(1.13) 

Here r > is a fixed macroscopic reference point and s G M with sT^/^ the 
longitudinal distance, s i— > A^is) is regarded as a stochastic process in s. 
Our central result is the convergence 

lim At{s) = kA{sk/2) (1.14) 

T — ^oo 

in the sense of convergence of finite dimensional distributions. The limit 
process A{s) is the stationary Airy process. Its scale is determined by the 
local curvature via k = ^26'4(r). The Airy process appeared first in the 
study of shape fluctuations for the polynuclear growth model jEj. It can also 
be obtained through edge scaling of /3 = 2, GUE Dyson's Brownian motion. 

2 Line ensemble, determinantal process 
2.1 Gradient lines 

In view of Figure it is natural to represent h in term of its level lines with 
the hope that they have a tractable statistics. In fact, it turns out to be more 
convenient to consider the gradient lines as drawn in Figure lU In Figure El 
the underlying lattice is distorted in such a way that the gradient lines be- 
come "trajectories" on a square lattice. It is this latter representation which 
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Figure 4: The gradient lines for the tiling of Figure\^ 



will be used in the sequel. Clearly, the surface statistics can be reconstructed 
from the statistics of the line ensemble. As first noticed by Okounkov and 
Reshetikhin ^2] the occupation number field corresponding to the line en- 
semble of Figure El has determinantal correlations. In this section we will 
rederive their results using the fermionic framework, which is a convenient 
starting point for our asymptotic analysis. 

The gradient lines of Figures EJ El are defined through 

t = j-i, he{t) = h{i,j)-e{i,j), (2.1) 

where 

i{i,j) = {t+j-\i-j\)/2 (2.2) 
labels the line, {i,j) G Z^. hi is increasing for t < and decreasing for t > 0, 

he{t) <hi{t + l), t<0, he{t) > he{t + l), t > 0, (2.3) 
with the asymptotic condition 

lim he{t) = -i. (2.4) 

t— >itoo 

By construction the gradient lines satisfy the non-crossing constraint 

h+i{t) <hiit~l), t<0, hi+i{t) <heit + l), t>0. (2.5) 



Height configurations h E are mapped one-to-one to gradient lines satisfy- 
ing (Q, (HI, and (123). 
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Figure 5: The gradient lines for the tiling of Figure{^ 



We extend he to piecewise constant functions on M such that the jumps 
are at midpoints, i.e. at points of Z + |. For a given hne, h^, let te^i < . . . < 
te,k{^) < be the left jump times with jump heights si^i, . . . , Si^k{e) and let 
< ti^k{i)+i < • • • < ti,k{£)+n{i) be the right jump times with jump heights 
—si^k(e.)+i, ■ ■ ■ , —se^k{e)+n{e)- It follows from ()1.4|) . ()1.5|1 that the weight for the 
line configuration {/i£}^=o,i,... is given by 



The line ensemble with weight fl2.6p can be thought of as world lines 
of "fermions", where t refers to time and j to space. It is then natural to 
introduce the random field of occupation variables, denoted by r]{j,t). Thus 
ri{j, t) = 1 if there is a hne passing at (j, t) and ri{j, t) = otherwise. As to be 
shown, the random field t] has determinantal correlation functions, which as 
one crucial ingredient relies on the non-crossing constraint ()2.5p . However, 
in previous applications only nearest neighbor jumps appear, whereas our 
model has the unusual feature that jumps of arbitrary size are allowed. 

2.2 Fermions 

The basic tool is the transfer matrix from ttot + l,tGZ. A fermion is 
created (resp. annihilated) at the position j G Z by the operator a* (resp. aj). 
The CAR algebra {a*,aj,j G Z} over Z is defined by the anticommutation 




(2.6) 



relations 



{ai,aj} = 0, {a*,a*} = 0, {ai,a*} = 6, 



(2.7) 



2 LINE ENSEMBLE, DETERMINANTAL PROCESS 



9 



for i,j G Z. First we consider t < —1, in which case only up-steps can occur. 
To each unit up-step at time t we assign the weight qt = which satisfy 

()2.6|) . The rule is that in a jump from i to j , j > i, one creates additional 
particles at sites m with i + 1 < m < j and annihilates particles at sites 
m with i < m < j — 1. E.g. if a fermionic world line jumps from —1 to 
3, one creates particles at positions 0, 1, 2, 3, and annihilates the particles at 
—1,0,1,2. This rule ensures the non-crossing constraint ()2.5p . since, if two 
fermionic world lines would intersect, a fermion is created twice at the same 
position, which leads to a zero contribution. The corresponding rule applies 
to t > with the difference that the jumps are downwards only. 
Let us define the operators 



(2.8) 



The transfer matrix from t to t + 1 is a sum of the n-step transitions T„ as 
f{t,t + 1) = 1 + qtTi + q^T2 + . . . + g^^n + • • • , (2.9) 

where 

(—1)" 

T„ = — -|— ttk, . . . ak„al^+i . . . al^^^. (2.10) 

The (—1)" prefactor results from the left ordering of the aj and a*'s. 

We would like to reexpress T„ in terms of products of the 6j's only. For 
n, m > the commutators are 

bnO'k ~ ^k^n "I" ^k+ni ^n^'k — (^k^n ~ (^k-n: [bny ^m] = 0, [b-n: ^-m] = 0. 

(2.11) 

These relations lead to 

di,...,d„>i j=i ^-^ ^ ^ 

di+2d2+...=n 

The Schur polynomials {pkiy)}k>o are polynomials such that 

exp ( ^t'^yfc j = ^Pi(y)t', y = yi,y2,..., (2.13) 

^ k>l ^ l>0 



and given explicitly by 



I Xj 

y 



Xl,...,Xl>l j = l 

Xi+2x2 + ...=l 
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Comparing with ()2.13|) yields 



f{t, t+l) = Y^ = exp ( J] qfA . (2.15) 

i>o ^ fc>l ^ 

We conclude that the transfer matrix is given by 

f(t,t + l) = exp(5^/l*+V^l|) (2.16) 
^ fc>i ^ 

for t G Z_ = { — 1, — 2, . . .}, and, by the same reasoning, 

f(t,t + l)=exp(5^g^(*+V^)^) (2.17) 

^ k>l ' 

for Z+. 

T{t, t + 1) is quadratic in the fermion operators. Hence it is the second 
quantization of a one-particle operator acting of £2- For easier reading second 
quantization is merely indicated by a "^"^ i.e. for A acting of £2 "we set 
A = T{A) as its second quantization. From ()2.16p . ()2.17|) we read off 



/ fc|m/2| \ 
T{t, t + 1) = exp ( ^ ^^"'^ ) 
^ fc>i ^ 



for t G Z_, and 

Mt+l/2) 



T(t, ^ + 1) = exp J] ^— ^ j (2.19) 



for t G Z_|_ with matrices defined through 



T{t,t + 1) are invertible with the £2-norms 
||T(t,t+l)|| < exp 



g|t+l/2| 
1 _ g|t+l/2| 

||T(t,t+l)-i|| < exp(^^A___j. (2.21) 

For the state at t = ±00 all sites in Z_ U {0} are filled, those in Z+ \ {0} 
are empty, which, together with the transfer matrices ()2.1(i|l . ()2.17j) deter- 
mines the Green's functions of an imaginary time (Euclidean) Fermi field. 
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It is inhomogeneous in space-time and uniquely given through the two-point 
function {a*{t)aj{t')) . To compute it correctly one has to employ the stan- 
dard finite volume approximation. We first restrict all world lines to lie in 
the spatial interval [— M, M\. Thereby the transfer matrix depends on M in 
the sense that all creation and annihilation operators with index \k\> M are 
set equal to zero. The state at ±00 is (1, . . . , 1, 0, . . . , 0)* which is 2M -|- 1 
long with the last 1 at site 0. The projector on this state is approximated 
through 

exp[/3iVM] (2.22) 

in the limit (3 ^ 00 with Nm = Y1\=-m ^l^k — Sfcli cil^k- We first compute 
the equal time Green's function through 

{a*{to)a,{to))T = (2.23) 

= lim lim lim Tr e'^^'^' n'^V/(t, t + l)a*a,- FT' fA/(t, t + 1) , 

M^ooL^oo/3^ooZ^,M,L ^ fto t l J 

where the trace is over the antisymmetric Fock space JF with one-particle 
space 12{\—M,...,M]). The products are time-ordered increasingly from 
right to left, which is indicated by the superscript t at the product symbol 
Y\. Zi3^M,L is the normalizing partition function, which is defined through the 
same trace with a*aj replaced by 1. As explained in Appendix \A.1\ ()2.23|) 
can be expressed in terms of one-particle operators as the limit M, L, (3 00 
of 



{a*{tQ)aj(to))T,i3ML 



'to-l L ^ 

1+1 \{'TM{t,t + l)e^''-'\{'TM{t^ + l) 



\t=-L t=to 

(2.24) 

Let P+ + P_ = 1 in £2 with P+ the projection onto Z+ \ {0} and let 
gGHght(to) ^ Y[ T{t, t+1), e^'=ft(*o) = Yl T{t, t+1), (2.25) 

t=to t= — oo 

and 

min(0,to)-l 00 

gGT(io)_ JJ T{t,t + 1), e^i(*°) = Yl T{t,t+1). (2.26) 

t=—oo t=max(0,to) 



By fl2.21|l the infinite products are well-defined, their inverses. The 
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T{t,t + l)'s commute and no time-ordering is required. Hence 
GT(to) 



mm(0,to)-l ^r|t+l/2| 



t=~oo r>l 



r>l 



with 



Y Y — ; — = Y''- (^o) 



t=max(0,to) 



„r/2„— rmin(0,to) „r/2 „r max(0,to) 



r(l — g^) ' ''^ r(l — g'') 

In ()2.24|) we take limits as indicated in ()2.2H|1 . Then 



(2.27) 



(2.28) 



{a* (to) a j (to)) T = [e^irft(*o)p_(P_e^"sht(to)gGi,ft{to)p_ + p^)-ip_eGHgM(to)j ^_ 

(2.29) 

(2.30) 



Let 



e^- = n Tit, t+1), = n ^(^, ^ + !)• 



t=— oo 



t=0 



Then e*^"ght(to)gGieft(to) = e^+e'^- = e'^-e^+ and, decomposing £2 = P-£2 
P+i2, we have 



a 
c 6 



.G4 



a' c' 
b' 



Thus 



P (P e'^^-'^sbtiU)) ^Gicftiio) p _|_ p^'j^^p 
and, since 



(aa')"^ 




fl-i 



-G+ 



(6')"' 



(2.31) 



(2.32) 



(2.33) 



we obtain 

g-G+p = P (^P gG'right(*o)gGlcft(<o)p _j_ P^^~^P 

Therefore 

{a*{to)aj{to))T = [e^irft(*o)g-G+p_g-G_gGHgM{to)]^^^ 
which rewrites as 

{aKto)a,{to))T = [e^T(*o)-G,(*o)p_g-(G,(io)-G,(io))] _ (2.36) 



(2.34) 



(2.35) 
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The Fermi field depends on T through q = 1 — 1/T. For this reason we keep 
the index T. 

Using the anticommutation relations ()2.7|1 in ()2.23|1 we immediately ob- 
tain 

(a,-(to)a*(to))T = [e^T(*o)-Gi(*o)p_^e-{GT(to)-G,(to))]^,^^ _ (2.37) 
Thus our final result for the equal time correlations reads 



l<0 

(«.(to)a*(to))T = ^[e«T(*o)-c,(*o)]^.^[,-G,(.oHG,(*o)]^^^. (2.3^ 



l>0 



To extend ()2.38p to unequal times we have to go through the same limit 
procedure as before. Since the argument is in essence unchanged, there is no 
need to repeat. We define the propagator from a to b, a < b, through 



b-l 

oG{a,b) _ 



JjT(t,t + l), e^("''^) = l, e^('''") = e-^('^'*). (2.39) 

t=a 

Using the identity 

e-^(°'*o)a^e«(°'*o) = b""^'''^] m.u^^ (2-40) 
for t >t', the full two-point function is given by 

l<0 

{aj{t)a*,it'))T = ^[e«T(o)-G,(o)+G(o,i)]^^^[g-GT(o)+G,(o)-G(o,0]^^^^_ (2.41) 



l>0 



2.3 Determinantal random field r]{j,t) 

Moments of the random field ri{j,t) introduced above can be expressed 
through fermionic correlations. We consider first equal time correlations. 
The basic identity is 

Et \flvTijk,t)] = (flal{t)a,,it)) , (2.42) 



\k=l / ^ k=l 
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where is the expectation with respect to the normahzed weight ()2.6|) . If 
{ji! • • • ! in} are distinct, then, as explained in Appendix IA.2t the fermionic 
expectation is determinantal and 



Et (^n^(jfc,^)j = Det(i?T(jfc,t;j;,t))i<fc,K„, (2.43) 

with 

RTit,t;j,t) = {a*it)a,it))T. (2.44) 

If coinciding arguments are admitted, then ()2.43|1 still holds with the con- 
vention 

^T(M;j,t) = ( , / r.wr.n (2.45) 

I {<{t)aj{t))T - = -(«iW<W)T for z > J. ^ ' 

()2.42|1 is easily extended to unequal time correlations. Let us consider 
n disjoint space-time points (ji,ti), . . . , {imtn) ordered increasingly as ti < 
^2 < • • • ^ ^n- Then the basic identity is 

Et {^Vijk.tk)^ = {a]Stn)a,Atn) ■ ■ ■a*^(ti)a,,(ti))T. (2.46) 
Using ()2.40|) the left hand side equals 

n 



ki,...,kn q=l 

ll,...,ln 

Let us set 

RT{j,t;j\t') = ( ; - (2.48) 

'-^ ^ [ -(ajv(i(:')a*(t))T. for t < t'. ^ ' 

Then the unequal time correlations are given by 

Et (j[vUk,tk?j = T^et{RT{jk.tk\]uti))i<k,i<n- (2-49) 

The identity ()2.49j) has been derived from left to right. One can read 
it also from right to left. Then Rf is the defining kernel, resp. Green's 
function, which is considered to be given and ()2.49|) defines the moments 
of some determinantal space-time random field over Z x Z. Of course, Rt 
cannot be chosen arbitrarily, since the right hand side of ()2.49|1 must be 
moments of a probability measure. For determinantal random fields over the 
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space coordinate only, compare with ()2.43|) . proper conditions on the defining 
kernel have been studied in detail [T7l[Tn]. The space-time variant is less well 
understood, see [Zj for a discussion. 

The determinantal property is preserved under limits. Thus through bulk 
and edge scaling further determinantal space-time random fields will be en- 
countered below. One of them is over Z x Z with equal-time given through 
the sine-kernel. The other is over Z x M with equal-time given through the 
Airy kernel. 



3 Limit shape 

On the macroscopic scale, in the limit T — »• oo, the random field i]{j,t) 
becomes deterministic with a profile given by 

r 1 forC<&-(r), 
p(C, t) = 1 ^ arccos (cosh(r/2) - e-^+l-l/V2) for 6^(r) < C < &oo(r), 
[ for C > &oo(r), 

(3.1) 

with 

b^ir) = -21n(l + b^ir) = -21n(l - e-l^l/^). (3.2) 

More precisely, for all continuous test functions / : ^ M with compact 
support 

T^H^^E/(^'/r,t/r)r/(j,t) = / dCdrp(C,r)/(C,r) (3.3) 

almost surely. ()3.3p assumes more spatial averaging than needed. In fact, it 
suffices to choose a test function whose support on the scale of the lattice 
diverges as T ^ oo and to properly normalize. 

As a consequence of ()3.3|1 the limit ()1.9j) holds, /ima can be read off from 
p.ip and is given in parametric form through 

h (uv)-l^ foT{u,v)eRl\V, 



where t = v — u and where C{u,v) is the unique solution ( of the equation 



2 



:u + v-\t\)= f (i_p(C',r))dC'-C (3.5) 



in the interval [^ool'^)' ^oo(t)]. While the limit ()3.3j) has been established by 
Okounkov and Reshetikhin compare also with Sectional the existence 
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of the limit shape has been proved before by Cerf and Kenyon .2j. Instead of 
(jl.Sp . they used the fixed volume constraint V{h) = 2C/j(3)T^, resp. V{h) < 
2Cij(3)T^, with the Riemann zeta function. They write the limit shape 
5*0, as a set of M.^, in the parametric representation 

So = {2(/(a, b, c) — In a, f{a, b, c) — In b, f{a, b, c) — In c)\ a, b,c> 0} (3.6) 

with 

-| /'2iT /'27r 

f{a,b,c) = — du dt;ln(a + 6e^" + ce™). (3.7) 
Jo Jo 

Here a, b, c denote the weights for the three orientations of the lozenges and 
/(a, b, c) is the corresponding free energy per unit area for the lozenge tiling 
of the plane. As expected from the equivalence of ensembles, the shapes 
given by ()3.4|) and ()3.6|) are identical. This can be seen as follows. Let 
z = {zi, Z2, Z3) represent a point on the limit shape. We compare Z2 — Z\ and 
2:3 — Zx (resp. z-i — Z2) for Z2 > Z\ (resp. z^ < Z\) for the parametrizations 
diH) and dini). This leads to a = 1, 6 = e'^/^ c = e'^/^ for Z2 > zi and to 
6 = 1, a = e"!"^!/^, c = e~^/^ for Z2 < zi. Since ()3.7j) is symmetric in a,b,c, 
one verifies that indeed 

/"^ (1 - piC, r)) dC = 2/(1, e-l-l/^ e-^/2) + C- (3.8) 

J-21n(l+e-l^l/2) 

According to ()3.4j) . /ima = on \ V. Close to the edge the height 
vanishes with the power 3/2. E.g. in the direction t = v — u one has 

/i^a(r, r) ~ cosh(r/4)7r-^2^/V/2 (3.9) 

with r the distance to the edge. The 3/2 power is known as Pokrovsky- 
Talapov law P^ . 

A limit shape theorem is a law of large numbers. It is available also 
for related tiling models. A famous case is the Aztec diamond 0. Cohn, 
Larsen and Propp |1] consider the 3D- Young diagrams constrained to the 
box aN X (3N x '^N with 7 ~ O {1) and compute the limit shape as 
00. In the line-ensemble representation their model corresponds to 
q = I with the boundary conditions that at t = —aN, j3N all lattices sites 
are occupied expect for those in the interval [1,7A^]. |lj computed the line 
density and from it the limit shape. Two or higher order point functions are 
not studied. From our representation we see that higher order correlation 
functions are determinantal even in this case. However the computation of 
the two-point function is more complicated, since one cannot rely any more 
on an expression like ()2.36j) . For a list of further limit shape theorems we 
refer to the survey [Til] . 
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The limit shape can be determined through minimizing the appropriate 
macroscopic free energy functional. The input is the microscopic surface 
tension at given slope Vh. For example in the (lll)-frame the surface tension 
c"(iii)(V/i) is given by ()3.7|1 . where a,b,c are defined through the prescribed 
surface tilt Vh. <J{iu) has been computed in [HI CHI H]- Correspondingly 
there is a surface tension in the (OOl)-frame, denoted by (T(ooi)(V/i). 

To obtain the free energy JF for some macroscopic height profile h over 
a bounded domain B, one argues that h is made up of little planar pieces, 
each one of them having the surface tension at the corresponding local slope. 
Adding up yields 



In our case we have B = M^, h is decreasing in both variables such that 
h{u,v) = OfoT {u,v)-^oo,a.ndV{h) = J^dudv h{u,v). The minimizer of JF, 
under these constraints and V{h) = 2C/j(3), is h^^, from ()3.4p . Equivalently 
one could minimize J-'{h) + V{h). 

Probabilistically, J^{h) + V{h) can be viewed as a large deviation func- 
tional in the sense that in the limit T — oo, with respect to the normalized 
probability Z^^q^^'^\ 



for given macroscopic height profile h j2]. 

Expanding 1)3.111) to quadratic order in 6h = h — /ima yields a heuristic 
formula for the covariance of the Gaussian shape fluctuations. In spirit it 
is proportional to (— — <9^)~^, hence like a massless Gaussian field. This 
implies in particular, that on the macroscopic scale shape fluctuations are 
small, of order InT only. Gaussian fluctuations are proved for the Aztec 
diamond in [Hj and for domino tilings of a Temperleyan polyomino in fTI] . 

The limit shape theorem ()3.1)1 implies that also the border step has a 
deterministic limit. We state a result, which is stronger than what could be 
deduced from ()3.1|) and which follows by the transfer matrix techniques to 
be explained in Section |S1 

Theorem 3.1. Let be the border step as defined in M.10\j . Then for any 
(5 > 0, c > 0, < < < oo one has 



lim P {\T-^bT{[uT]) - b^{u)\ > cT-^/^+^,u. <u<u+) =0. (3.12) 




(3.10) 




) (3.11) 
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4 Bulk scaling, local equilibrium 

For local equilibrium we zoom to a point (Co, to)T with 6^(to) < Co < &oo(to) 
at average density p — p(Co, To), which means to consider the random field 

V^'^'^j, t) = viiCoT] + j, [toT] + t) (4.1) 

with {j,t) e and [ ] denoting the integer part. Properly speaking we 
should keep the reference point (Co, Tq) in our notation. Since it is fixed 
throughout, we suppress it for simplicity. In the limit T oo, Vr'^^U^t) 
becomes stationary. Then at fixed t, one has to fill the Fermi states up to 
the density p which implies that r]^^^{j,t), t fixed, is a determinantal point 
process on Z as defined through the discrete sine-kernel. Only at Tq = 0, the 
inhomogeneity of the underlying 77-field can still be seen, which, of course, is 
an artifact of our coordinate system. In the (111) projection the line Tq — 
would be just like any other local slope with a corresponding stationary 
distribution of lozenges. The case tq = can also be treated. For simplicity 
we omit it and require To 7^ 0. 

Let us define the kernel S{j,t] j' ,t') by 

S{3, t; f, t') = ^g^ii^ / dA; exp [zk{j' - 3) + {f - t) ln(l - e-\^^\/\-^')] 

(4.2) 

for To > and 

S{3.t-j\t')^^-^^^^ I dfcexp [^k{f -3)-{t' -t)Hl-e-\^^\l\^^)\ 

(4.3) 

for To < 0, where 

^^'^^ " \ [7rp,27r-7rp], if t < f , 

and sgn(t — t') = 1, if t > t', and sgn{t — t') — —1, ii t < t'. In particular at 
equal times 

= (4.4) 

which is the sine-kernel. S depends on the reference point (Co, To)- In the 
particular case of equal times the dependence is only through the local den- 
sity. 

Theorem 4.1. In the sense of convergence of local distributions we have 

hm VT'\j,t)=v''''%J,t). (4.5) 

1 — >oo 



4 BULK SCALING, LOCAL EQUILIBRIUM 



19 



Fottq > 0, r]^^^'^{j,t) is the determinantal space-time random field with defin- 
ing kernel and for tq < the one with the kernel d^.^p . 

Remark: Theorem 14.11 is identical to Theorem 2 of \l2i- We use here an 
integral representation for the defining kernel Rt which differs somewhat 
from the one of ^2] and which turns out to be convenient in the context of 
the edge scaling. 

Proof: We consider the case tq > only, since tq < follows by symmetry. 
Let us set 

BtU, t; f, t') = e^(^-)-^(^-')i?T([CoT] + j, [tqT] + t; [CoT] + j', [tqT] + t'), (4.6) 

where Rt is defined in ()2.48j) and g{j) = j |ro| Tln(l — l/T)/2. The determi- 
nant in ()2.49|1 does not change under similarity transformation, in particular 
not under multiplying by e^^"*^"^*^"^^. Therefore 

Et ^VT'\jk,h?j = Det {BTUk,h;jiMi<k,i<m (4-7) 

and we need to prove that pointwise 

hm BTU,t-f,t') = SU,t;f,t'). (4.8) 

T — >oo 

First consider t >t' . For tq > we take T large enough so that tqT + t' > 
(this simplifies ()4.1Hj) below). Using ()2.49j) we obtain 



BTU,t;f,t') = e3(^>5(^') (4.9) 



X ^ [e'^T(^oT)-Gj,(ror+t)j^ ^ ^^^Gx(roT+t')-Gt(roT)j^ 



KO 



An explicit expression for the matrix elements of the two-point functions can 
be found using the translation invariance of the one-particle operators. In 
Fourier representation they are given by 

exp ( aj.arj ~ 2~ I ^■^P(~^^(^ ~ "^)) f cTrC*'^'^ j dk 

(4.10) 

for Ot G M. Then using ()4.1Up and changing / into — /, we have 

BT{j,t;j',t') = ake''(k)T^^,(k,t)^-ik(CoT+j)^-iki 

27r ] 
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where 

a{k) = il-q)J2 r) (^''' - r^'^'-'^e-^'^) (4.12) 

r>l ^ ' 

and 

To study the asymptotic of integrals as 1)4.1111 we consider the complex k 
plane and regard the integration in ()4.1H) as being along the real line from — vr 
to TT. Such a line integral can be deformed to another path C with the same 
endpoints. The complex integration along C will be denoted by dfc ■ ■ ■ . 
In the particular case when the path is on the real line, say from a to 6, the 
integral will be denoted by j^dk ■ ■ ■ . 

Let us consider the following four paths: = — — >• vr, ^i(p) = —n+ip 
7i + ip, ^2 = — TT — > —TT + ip, and ^3 = n + ip n with < p < Tq. The factors 
in ()4.11|) are integrals along ^q. Their integration contour can be deformed 
from ^0 to ^2 without changing the integrals, since the integrands are 

holomorphic. Moreover the integrals on ^2 and ^3 cancel exactly because of 
periodicity of the integrands. We transform the integral in k into the integral 
along .^1(6') and the one in k' into the integral along ^i{6 + e), with < e -C 1 
and 6 = — roTln(l — l/T)/2. 6 is chosen such that the exponentially large 
function in T passes through the critical point of a{k'). Consequently we 
have 

g(j) — Qij'^ r r 

BT{j,t-f,t') = \ / dk' / dke'•^i>'^-''^>''))T^^CoTik'-k) 

X g'/'9('=.*)-'/'9('='.*')g*(fc'i'-fci)(']_ _ g^Cs'-fe)")-!^ (4.14) 

As T 00 we obtain 

-rro/2 

(t(A;) = 2i — sin((A; - iro/2)r) + O {l/T) . (4.15) 

r>l ^ 

Therefore the terms that increase or decrease exponentially in T in fl4.14|) 

are E{k) and —E{k'), where 

E{k) = 2i 2^ — sm{{k - iro/2)r) - iCok. (4.16) 

The critical points of E{k) are 

/ g-Co+ro/2\ 

±kc + iTo/2, kc = arccos ( cosh(ro/2) j e M. (4.17) 
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Figure 6: Deformation of integration paths. The original integral, along Iq, 
is deformed to the sum of integrals along Ii, I2, and J3. k is integrated along 
the dashed lines and k' along the solid lines. The full dots are the critical 
points of a{k). 



For Im(fc) = ro/2, Re(£^(/c)) = ^0^0/2, the analysis of Re(£'(fc)) for k close 
to the line Im(/c) = ro/2 shows that, for Re(i?(fc)) G [— vr, — /cj U [fcc,vr], it 
decreases when increasing Im(A;) and, for Re(i?(A;)) G [— fcc, /cc], it decreases 
when decreasing \m.{k). 

Next we transform the integral into a sum of three terms, the first two 
vanish as T ^ 00 and the third one gives the final result, see Figure IHl We 
have /^^ dk' dk - ■ ■ = f^ ^ dk' d k ■ ■ ■ + fj^ dk' dk ■ ■ ■ + Jj^ dk' dk ■ ■ ■ , where the 
integrand is the one of fj4.14|l . Let us compute the three integrals separately. 
For the integration along I3, first we integrate out k taking the residuum at 
k = k'. Then changing the variable to z = k' — i6 we obtain 

[ dk'dk--- = — / ' (ize'^ii^+'^,t)-v^,(^+ie,t')ei<f-^). (4.I8) 

The asymptotic of ipg{z + i9, t) is 

ip{z, t) = lim ifig{z + i9, t) = -t ln(l - e-^^/^e-''). (4.19) 
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The integrals along Ji and I2 are treated in the same way. Let us estimate, 
e.g., the one along Ii. First we integrate in k. The integral | J dfc ■ ■ ■ | , 
such that the integration avoids the two arcs of circle of radius e around 
the critical points (see Figure El), is bounded by O [e-"^^ /{eT)) for some 
a > 0. O (e~"^"^/T) comes from integrating e^^''^^ and O (l/e) because the 
minimum of \k — k'\ equals e. The integration through the two arcs around 
the critical points is bounded by O [e"" / (eT)) for some a' > 0, because 
the integrand is at most of O (e"'^"^/e) for some a' > and the length of 
the path of integration is 0{1/T). We choose therefore e = l/T, so that 
I J dfc ■ ■ ■ I < O (1). The integration in k' gives an extra-factor (9 (l/T), and 

lim / dA;' dA; ■ ■ ■ = 0. (4.20) 
Summarizing for t > t', we have proved that, 

lim BT{j,t;f,t') = — dze^("'*-*')e^"(j'-^), (4.21) 

where p = kc/n and ip{z,t) as given in ()4.19|) . In particular for t = t', 
(p{z,t — t') = 0, which implies ()4.4|) . The case t < t' is treated in a similar 
way, leading to 

Mm BT{3,t;j\t') = / dze^^'''-''^e"^^'-^\ (4.22) 

27r 

Therefore 



\im^ET[l[vT'''Hjk,h)j = Det(5(jfc,tfc;j,,tz)),<,, 

= Ebfnr/^''^U,4)). (4.23) 



The proof for tq < is identical. □ 

()4.21|1 . ()4.22j) define a space-time homogeneous Fermi field. Physically it 
corresponds to fermions on the lattice Z in their ground state at density p 
and with kinetic energy 

Ekin(A;) = ln(l - e-l"ol/2-^'=«s>^-o)_ (4 24) 



i^kin is complex reflecting that the fermions have a drift. 
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The moments ()4.23|) define a probability measure Pb on the lozenge tilings 
of the plane, where the relative fraction of their type depends on the reference 
point (Co, To)- Pb is a Gibbs measure in the sense that its conditional expec- 
tations satisfy the DLR equations. We refer to [Sj of how DLR equations are 
adjusted in the context of surface models. Pb is translation invariant with a 
definite fraction of each type of lozenges. Pb is even spatially mixing, since 
truncated correlations decay to zero. One would expect that Pb is the unique 
Gibbs measure with these properties. A proof would require that the same 
limit measure Pb is obtained when other boundary conditions are imposed, 
at fixed lozenge chemical potentials. To our knowledge, only for the surface 
model studied in |3] such a uniqueness property has been established. 



5 Edge scaling 

For the edge scaling one zooms at a macroscopic point lying exactly on the 
border of the facet, i.e. at {(o,tq)T with (q = 6oo(to)- For simplicity we set 
To > 0. To < follows by symmetry. Since at the edge the step density is 
zero, one has to consider a scale coarser than the one for the bulk scaling 
in Section |3] From our study of the PNG droplet we know already that the 
longitudinal scale is T^/^ and the transversal scale is T^^^. On that scale the 
curvature of 6oo cannot be neglected. Therefore the correct reference points 
are 

t{s) = [toT + sT'/% (5.1) 
j{r,s) = [6^(ro)T + 6'^(ro)sT2/3 + X(ro)3^TV3 + rTV3]. 

Note that (r, s) G M?. The discrete lattice disappears under edge scaling. 
Let us abbreviate 

«i = &oo(ro) = -21n(l-e--o/2), 

«2 = -b'^iro) = e-W7(i_e-<'/2), (5.2) 

«3 = bliro) = e-W72(l-e-W2)2. 

Then the edge-scaled random field reads 

rj'^^^%r,s) = T^/^r]T{[a,T-a2sT'/^+'^ass'T'/^+rT'/% [toT+sT^/^]). (5.3) 

The prefactor T^/^ jg ^]-^g volume element for rT^^^. Properly speaking we 
should keep the reference time tq. Since it is fixed throughout, we suppress 
it in our notation. 

Since rj^^'^ is determinantal, so must be its limit. For the PNG droplet 
under edge scaling the limit is the Airy random field and, by universality, in 
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our model the steps close to the facet edge should have the same statistics in 
the limit T ^ oo. The Airy field is determinantal in space-time with Green's 
function 

K^^^{r, s; r', s') = sgn(s' - s) [ d\9{\{s - s'))e^^''-'^ Ai (r - A) Ai (r' - A) , 

(5.4) 

where the step function 9{s) = 0, if s < 0, and 9{s) = 1, if s > 0. The Airy 
field is stationary in time. In particular, the equal time correlations are given 
through the Airy kernel 

^0 

K^"y{r,s;r',s) = / dA Ai (r - A) Ai (r' - A) (5.5) 

J — oo 

= — ^(Ai(r)Ai'(r ) - Ai(r')Ai'(r)). 

Theorem 5.1. Under edge scaling i5.^) the correlation functions have the 
following limit, 

uniformly for r^, Sk in a hounded set. Here k = ^26'4(to). In particular for 
the process ri^^'^{f,s) = J dxf{x)ri'^^'^{x,s), smeared over continuous test 
functions / : M — > R with compact support, one has 



lim Vrif^s) = / dxf{Kx)r]^''^x,SK/2) (5.7) 

in the sense of the convergence of joint finite-dimensional distributions. 

To prove Theorem 15.11 one only has to establish that under edge scaling 
()2.48|1 converges to ()5.4j) . We define the rescaled kernel ()2.48j) as 

KT{r,s-r',s') = ^-^^T'/'R^{j{r, s),t{sy, j{r' , s'),t{s')) (5.8) 

where c/(r,s) = -j(r, s)(roTln(l - l/T)/2 + sT'^/^ln{l - l/T)/2) and 
RT{j,t;f,t') from 

Proposition 5.2. The edge-scaled kernel i5.^} converges to the Airy kernel 
hm Krir, s; r' , s') = k'^K"^^^ f-, ^s- ^s'] (5.9) 
uniformly for r, r', s, s' in bounded sets. 
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Granted Proposition 15.21 we establish Theorem 15.11 
Proof of Theorem I5.lt From ()2.49|1 and ()5.Hp it follows that 

Et (j[vT''\n,Sk)^ =Det{T'/'RTU{n,Sk),t{sk);j{n,si),t{si))) 

(5.10) 

This determinant does not change when multiplied by the factor e ,s ) 

and therefore 



l<kl<m ■ 



Et |^n'7?''('"fc,Sfe)j = Bet {Krirk, Sk;ri, si))^^i^ i^^. (5.11) 
Note that g{r, s) diverges as T — oo. On the other hand 

, (n (^, f ..-))) - (^f . f ■ 

(5.12) 

Theorem 15. II thus follows from ()5.9|) . □ 

We turn to the proof of Proposition l5.2l As bounded set we fix throughout 
a centered box B C M^, where the dimension d depends on the context. 

Proof of Proposition 15. 2t Let us first consider S2 > si. By definition 
of KT{r2, S2] ri, si), ((2311), (F^ - and (jQEI), we have 

g-3(n,si) 

X /''dA;e-*^^'('^^'^^)e*'='e^">i('^"^*"-'^"^""")eS">i^"^""" (5.13) 

1 /"TI" , , , 



27r 



where /i„ = q"'/'^ /n{l — q"-), i/„ = iin<f"^°^ , and v?^ = z/.„(l — g"'*'^^''^). As in 
Section |3] we regard the integrals in ()5.13|) as complex line integrals and use 
the notation explained below ()4.13j) . 

The integrands in ()5.13|1 are holomorphic away from {A; G C | 'Re{k) = 
0, I lm{k—iTo/2) \ > Tq/2} and the straight path from — vr to vr can be deformed 
provided no singularities are touched. In our choice the deformed path has 
three straight lines, the first one from — tt to — tt + ij3i(T), the second one 
from —71 + iPi{T) to vr + iPi{T), and the last one from n + iPi{T) to tt with 
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Im 

^0 




Figure 7: Deformation of the integration path. The original integral, from 
—TT to TT, is deformed along the integral on the dashed path. 



Pi G (0, To), see Figure To be precise, the path along the real line touches 
at k = the starting point of a branch cut of the term in the exponential, 
but still the integral remains unchanged by the above deformation. Since the 
integrands are 27r-periodic along the real axis, the first and the last integrals 
cancels exactly. l3i{T) is determined such that the terms in the exponential 
are purely imaginary. We obtain 

Am = -\ (s.T2/3ln(l - 1/T) + roTln(l - 1/T)) , z = 1,2. (5.14) 

We also define / = LT^/^. Then the summation goes over L E T^^/^(Z_U{0}) 
and 



p-9iri,si) 2^1/3 _ _ ^ ^ 



where 

and 
ML) 

ML) 



g-g(r2,S2) 47;-2 

Ler-i/3(z_u{0}) 
6i=j{r„s,)f3i{T)-f3,{T)LT'/^ 



2i /i„ sin(A;n)e 



(T)n 



n>l 



(5.15) 



(5.16) 



dk, (5.17) 



g.fc',(r-2,S2)e-.fc'LTV3^^p 



2i fin sin(/c'n)e~ 



n>l 



dk'. 
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Finally defining Ji{L) = T^^^Ji(L), we have 

LeT-i/3(z_U{0}) 

(5.18) 

For the case S2 < si the result is 

Kr{r„s,■,n,s^) = - (47r^TV3)-ie^^(---)(i+^(^'^)) Ji(L) J2(i.). 

LeT-i/3(z+\{0}) 

(5.19) 

Now we proceed as follows. First we prove that, as T — > oo, 
Ji{L) — > ^ Ai {^-^^) for L G i3, by using the steepest descend curve for the 
term which is exponentially small in T. Secondly we consider separately 
S2 < Si and S2 > Si. In the latter case, for large L, we need the steepest 
descend curve for the whole integrand. The same strategy has been used 
in p. In the case S2 < si, for large L, the steepest descend curve does 
not exist anymore. On the other hand the term e"^'-'^^"''^-' serves as a 
convergence factor and we only need to find bounds for the Ji{L). 

Convergence for L in a bounded set. 

Let L ^ B. The integral Ji{L) is written as 

Ji(L) = T'/' r e^'^i.-We^'^^^''' dfc, (5.20) 

J —IT 

where 

^i^y(A;) = -tkT^^jin, Si) + 2tY ^e"^'" sin(A;n). (5.21) 

We make a saddle point approximation by using a curve which, for small k, 
is very close to the steepest descend curve for ipi^k), where 

^p{k)= lim {ikLT^/'^ + ^i^T{k))/T (5.22) 

T — ^oo 

and the convergence is uniform for (si,ri,L) G B. For the limit we obtain 

i){k) = + 2ik In (l - e"^"/^) , (5.23) 

where 

2isin(fcn) —. (5.24) 

n>l 

In particular tp{k) is holomorphic in C \ {k = x + iy & C \ x = 0, \y\ > tq/2} 
and the whole integrand is 27r-periodic along the real axis. 
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Figure 8: Deformation of the integration path. The path from —7i to tt is 
deformed into C plus the dashed ones. 



Instead of integrating along the straight path — tt ^ tt we integrate along 
C = {k = x+iy, y = — \x\ / VS}, see Figure|Hl For x small this path is almost 
at steepest descend. The real part of il>{k) reaches its maximum at = 0. 
To evaluate the errors for x away from zero we prove that the real part of 
ip{k) is strictly decreasing for \x\ increasing. By symmetry we consider only 
X E [0,7r]. A simple computation gives 



d^(fc) 
dx 



+ 1/^3) ln(Q) 



with 



and 



Q 

ReQ 
ImQ 



'I _ e--o/2)2 



(5.25) 



(5.26) 



cosh(ro/2) — cosh(x/ y/S) cos(x) 
2sinh2(ro/4) 
sin(x) sinh(x/ v^) 



2sinh^(ro/4) 



(5.27) 



Using that cosh(x/v^) cos(x) < 1 and is maximal at x = 0, we have 
ReQ{x) > Re (5(0) = 1, the inequality being strict if x 7^ 0. Obviously 
ImQ < 0. Therefore 



1 



2V3 



ln((Reg)2 + (Img)2) + arctan(Img/Reg) < 

(5.28) 
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for all X G [0, tt] and for all tq G (0, oo). The inequality is strict if a; 7^ 0. 
Since ^MM = Re (^^) and by ReV(fc) is maximal at = 0, 

V'(O) = 0, and is strictly decreasing for increasing. 

Let us fix e, < £ ^ 1, and let be the part of C with x G 
Then the contribution at Ji{L) coming from C \ is exponentially small in 
T. 

Lemma 5.3. For some 5 > 0, 

ML) = 0{e-'^)+T'/' [ e'^i'-W^e*'^^^''' dA;. (5.29) 

Proof: Let be the part of C with x G [e, vr] and C~ the one with x G 
[-TT, -e]. For x>e, Reipik) < ReipiO) - 25 < for suitable 6 = 5{e) > 0. 
In addition 

^i,T{k)T + tkLT^/'"^ = ^{k)T + 0{{L- n)T^/^ + siT^/^) . (5.30) 

Then 

^t^^.T{k)T ^ikLT^ri ^ ^^ST r 2 ^(Re^(k)-S)T^O(iL-n)T^/3+sjTy3) 

~ Je VS 

(5.31) 

For (L, Si,ri) G B, the integral on the right side is uniformly bounded and 
therefore 

leM^^^e'^'^T^^-Uk =0{e-'"'). (5.32) 

Similarly for the integral along C~ . □ 
Lemma 5.4. Uniformly for (L,ri,Si) G one /ias 

Ji(L) = O (e-^^) + O (T-V3) + ^ Ai f ^i— (5.33) 

K \ K J 

for large T , with k, = 

Proof: By Lemma f5. 31 we have to evaluate the contribution of the integral 
along Ce. For k close to we have 

iJi,T{k)T + ikLT^/^ = -\ia^k^T -ikT^'^{rx- L) (5.34) 

+ O {sjk + siPt2/3 + k^T) . 
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Let be the part of C with x G [0,£:] and the one with x G [— £,0]. 
Then ■ ■ ■ = ■ ■ ■ + ■ ■ ■. We consider exphcitly only one of the two 
integrals, the second being evaluated in the same fashion, 

2^1/3 f ^^i,T(k)T+ikLTy^ 





e 


§jQ3 








L 


e~ 











= / e-i^"^'='^e-^'=^'^'('-^-^)dA; + Ei(L). 

Jc+ 

The error term is the integral along C+ with integrand 

2.1/3g-|m3fc^Tg-ifcTi/3(r.i_L)^^0(s2fc+sifc3T2/3+fc5T) _ ^-^ (5.36) 
= 2^1/3g-|m3fc='Tg-ifcTi/3(,.i_L)gCi(s2fc+.ifc3r2/3+fc5j.)^ ^^2^ _^ S-^PT^/S + Pt). 

The term in the exponential is — "^ia^k^Ti^l + xi) — ikT^^^{ri — L)(l + X2), 
where xi ^i-nd X2 can be made arbitrarily small by taking e small enough (si 
is bounded). With the change of variable z = kT^^^ we obtain 

EJL) = I e-*l"3(l+Xi)^^-i(l+X2){ri-L).(-) fg2 ^ ^^^3 ^ ^5y-l/3^ 

(5.37) 

Remark that at the boundary of the integration, the real part of the integrand 
behaves as e~3"3e t ^jj-^jg integral is uniformly bounded in T for (L, ri, Si) e 
B. The same holds for the integral on C~. Consequently Ei{L) = O (T~^/^) . 

Next we extend the integration from to — 7rT^/^(l, cos(7r/6)) and 
7rT^/'^(l, — cos('7r/6)), obtaining the path Di. In this way we add an error 
of O (e^'^''-^-'"^) with 5'{e) ~ e^. Similarly we can complete the path up to 
X = ±NttT^^^, y = —N7iT^^^/\/3 by straight lines. The integral is equal to 
the integration from —NttT^^^ to NttT^^^, since the function is peri- 
odic in the real direction and the error added by completing the integral is 
exponentially small in T, for all N. Therefore we may take the limit N 00. 

Finally we obtain ()5.33p . since 



re-^§-^^-(--^)dz = -Ai ('-1^] 

J-00 ^ V ^ / 



(5.38) 



with K = \/2a^. □ 
Convergence of K'r{r2, S2',ri, si) with S2 < Si. 



5 EDGE SCALING 



31 



Lemma 5.5. Uniformly for (si,rj) G i3, z = 1,2, 



(5.39) 



with K = v2a3. 



Proof: Since (ri,r2) G i3, let us set Lq such that Lq < 2{\ri\ + \r2\ + 1) for 
all ri, r2. -R't can be transformed into an integral adding an error 
Let us fix an e with < e <^ 1. Then 



r ° „-7.A ^-ej-'- e"^'^ 

-i^T(r2,S2;ri,si) = / ML)J2iL)-—dL+ / Ji(L) J2(L)— ^ dL 

Jo 47r^ Jig 47r^ 

/•oo p-LX 

+ / J,(L)J2(L)^dL + 0(T-i/3) (5.40) 

with X = i(si - S2)(l + C(l/T)) > 0. Since \ Ji{L)\ < T^''\ i = 1,2, the 
third term is bounded by T^/^e~^^^''^^/X — as T — oo. By Lemma 15.41 
the first term converges, uniformly for (wj, Sj) G i3, to 

^° . /ri - L\ . /r2 - L\ .£2^1 







as T — >■ oo. 

We consider the second term. We have already established the point- 
wise convergence of Ji{L) to ^Ai(^^^^). If we obtain that for large T, 
\Ji{L)\ < G with a constant G independent of r^, and L G [Lq, eT^/^], then 
by dominated convergence 

(5.42) 

uniformly for (r^, Sj) G B. This property is proven in the following lemma. □ 

Lemma 5.6. For L G [Lq^eT'^/^], \ Ji{L)\ < G with the constant G indepen- 
dent of Si,ri, and L, provided < £ ^ 1 and T large enough. 

Proof: The exponential terms in ()5.20|1 are purely imaginary for k real. Let 

us set 

i,{{k) = -{ikLT-^l^ + ^i,T{k)), (5.43) 



then 



J — 7r 
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In particular for k close to 0, 



^Pl (k) = -^a^e {l + O {e + s,T-'/')) - k{n - L)T-'/' {l + O (s?T-V3)) 

(5.45) 

Since {n - L)T-'^'^ ~ 0(e) at most, we set L = {L - ri)T-'^'^. i/j{{k) has 
two local extrema at ±k{L) with 



(5.46) 



and Co = (2^3) Moreover for \k\ > 2k{L), il){{k;) is strictly decreasing. 

Ji(L) = ■ • ■ = Eti 4 ■ ■ ■ where A = [-vr, -2c^Vl], h = [-2c^Vl, 0], 

J3 = [0,2co\/i'], and 1^ = [2co\/l, vr]. The integrals along Ji and /4 are 
evaluated similarly and so are the integrals along I2 and J3. We present in 

detail only the integration along J3 and I4. Let 7 = ^/l. Then 



/4 



2^1/3 / e*'^i^'')'^dfc = T^/^ 

^2co7 



?i(7r) 



/(M)e™' dw 



(5.47) 



m{2co7) 



where u = ipi{k) and /(m) = ^''^^^ ■ Integrating by parts we obtain 



rp-1/3 



f{n)- 



iuT 



iT 



u{n) 



u(2co7) "'«(2co7) 



u{tt) 



iuT 



df{u) e 



du iT 



du. 



(5.48) 



For k e h with \k\ < e follows from that ^ < and ^ > 0. For 

k > e, 



du ~ 



1 + 6-""° - 2e-"""/2 cos(A;) 



1 + e-^o - 2e-^o/2 

Then for k e h with A; > £, 4^ < and > 0. Therefore ^ 

du 



(5.49) 



-(^)-3g| where ^ < and g| > for every point in h. Thus ^ does 
not change sign along J4, and 



-f4 



<^(l/(«W)l + l/W2co7))l). 



(5.50) 



Using ()5.49|) . for T sufficiently large 



2 ln(l - e-""/') - 2 ln(l + e-^^'^) +7^ + (siT-^/^) | ' 



< I ln(l - e"^"/^) - ln(l + e"^°/^) 



-1 



(5.51) 
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provided e small enough (which implies 7 sufficiently small). The second 
term is bounded by 



l/(^(2co7))l 



1 + C» (72 + siT-Vs) 



< 2/7^ (5.52) 



for e small and Si e B. Therefore we have, uniformly in (ri, Si) e B, 









< 


J 1 4, 









2Gi 



2 2Gi 2 



T2/3 ^ (L - n) - T2/3 ^ (Lo - n) ■ 



(5.53) 



Next we estimate 



r /-Scot /•C27 

/ ■ • • = / e^^i'^'^)^ dA; = T'/' / e^^^'^)^ dA; (5.54) 

Jig Jo J —ci"f 

where ^jj{k) = ipKk - k{L)), a = Co(l + 0(7)) and C2 = Co(l + 0(7)). Let us 
define the paths ^0 = {k = x,x : — C17 C27}, ^1 = {k = — Ci7e~*''', ip : ^ 
/4}, ^2 = {k = e— : -C17 ^ C27}, 6 = = C27e*'^,95 : 7r/4 ^ 0}. 



Then /^^ • • • = /^^ ■ ■ ■ = Xli=i 4 • • • ■ The integrals along ^1 and ^3 are 
estimated in the same way. 



and therefore 

J^i Jo 



dip. (5.56) 



Since '0(A;(</?)) = -0(0) + |<(0)A;((/?)2(1 + 5i((^)) with Si{cp) ^ as £ ^ and 
k{ipY = 0^72(1 + O (7))e-2^v^ one has 

ImV^(^((^)) = -i<(0)(^((^))2(l + 52(<^))c^f (l + (^(7))sm(2(^) (5-57) 

with S2{(p) — > as £ — > 0. Moreover, for e small enough, sin(2(/7)(l+52((/?))(l+ 
C (7)) > Prom this it follows 



■7r/4 



< 2rV3co7 / e^'=o7^^"(o)W2d<^ (5.58) 
Jo 



TC272 ^"(0) 
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We compute <(0) = {l + O (7^ + s^T-^/^)). Therefore for si E B 

and T large enough, 



< 



< 



4 



(L-n) - {Lo-nY 

Next we need to evaluate the integral along ^2, 



Then for sufficiently small S, 



C27 



-C17 



tV3 / e^V^WT^^ 



/C27 
-C17 



<j.l/3 / grV"(0)x2(l+<5)/2^^ 



6 

00 J — t 



-00 



Thus we have, uniformly for si e B and T large enough. 



/ 



< 



8 



+ 



(i^o - ri) ^Lo-ri' 



Therefore Ji{L) is bounded by 

\ML)\<^ + 



20 



+ 



20F 



(5.59) 



gr^"(o).V2(i+oW)d^_ (5.60) 



(5.61) 



(5.62) 



(5.63) 



Since Lq — > 2 and X > 0, it follows that Ji(L) J2(L)e is bounded by 
an integrable function on [Lq, eT^/^] for < £ <S 1 and T large enough. □ 

Convergence of KT{r2, S2;ri, si) with S2 > si- 
Lemma 5.7. Uniformly for (sj, r^) e i = 1, 2, 



lim i^T(r2,S2;ri,Si) 

T— +00 



ri — L 



K 



Ai 



r2 — dL 
2" 



(5.64) 



with K = v^2q;3. 
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Proof: Let us set Lq such that Lq > 2(|ri| + \r2\ + 1) for all (ri,r2) G 
Then the sum in Kj- can be approximated by an integral at the expense of 
an error O (T^^^^Y Let us fix e, < £ ^ 1. Then 



KT{r2,S2;ri,si) 



-,LX 



-Lo 



J,{L)J2{L)^dL + 



-Lo LX 

ML)J2iL)-^dL 



+ 



-,LX 



Ji{L)J2{L)^dL + O , (5.65) 



with X = i(s2 — Si){l + O (1/T)) > 0. The convergence of the first term 
has already been proved. Let us set L = — (L — ri)T~^/^. In the remainder 
of the proof we set 

^ ~ (5.66) 



^P^k) = -^piTik) - kL. 



First consider L > e. 



(5.67) 



With the change of variable u = ip{k), f{u 
we have 



^^^i^, and integration by parts, 



max 

T k<=[ip{-TT),ip{Tr) 



dfiu) 



du 



To compute 



df{u) 
du 



we use 



d/(n) 
du 



' du\~3 d^u 
^ dk I dfc2 • 



du 
dfc 



(5.68) 
is flOni) with 

L replaced by —L. It is easy then to see that uniformly for si G 
maxfce[^(_^),^(^)] < GiL~^ for Gi = 2/(sinh(ro/2)£)2 < oo. Then 

for a suitable constant G2 < oo, 

\ML)\<G2{n-Ly\ 



(5.69) 



The same holds for J2, therefore the third term in ()5.65|1 is bounded by 



Gl 



dL, 



(L+ Inl + |r2|)2 

which is convergent for T finite and vanishes for T ^ 00. 
Finally we consider < L < e. Let us set P 

integrate over C 



(5.70) 



^2(cosh(ro/2) - 1). We 



k = X + iy{x),y{x) = —\/y{Oy + x^/Sj, with iy{0) the 
stationary point of tp{-,L), see Figure El y{0) = —i3\/l + O [L^^^j and C 
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Im 








To 

2 




— TT 







TT 


1 


TO 

2 


^ 

Re 



Figure 9: Deformation of the integration path. The path from —7i to tt is 
deformed into C plus the dashed ones. 



is almost the steepest descend curve for x small. This path has the property 
that the real part of iplk) is strictly decreasing as increases and 



^Pi^ym = -l^L'^' (l + O (s?T-V3 + L 



(5.7i: 



We divide the integral in the part with |x| < e and the remainder, 

J--K Jc JCe 

(5.72) 

where 



ST-mr,-Lf/^ 



Jc\Ce ^ 

(5.73) 

We then need to integrate only close to x = 0. We first establish some 
properties of tpik) for x = 0. 



Lemma 5.8. 

■■\ dip{k,L)\ 

dfc _\k=iy(&) 

.••\ d^4)(k,L) I 

^^^) dfc2 \k=iyiO) 



/3 



-Wl + o(l + VlsiT~'/^], 



(5.74) 



^t + O L + s^T-'/'). 



Proof: i) follows from Equation 1)5.711) and ii) because k = iy{0) is a sta- 
tionary point of iIj{-,L). iv) follows from ()5.34)1 because 2^3 = Finally 
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let A = V L. Then 



(5.75) 



dA dfc2 dfc3 dA 

and evaluating at = iy{0) and A = we obtain iii). □ 

With these properties 
Ji(L) - E2{L) = T^'^ [ e^^'^)^ dk (5.76) 



e 



o(L{k~iy{(d)fT+\/l.{k~iy{0)fsiT^/^+L{k-iy{0)fT+{k-iy{^^^ 

xe V / . 

Let 7 = v^ri — L, then \/l = 7T~^/^. Let A;' = k—iy{0), then the integration 
is along C'^ = Ce + iy{0). 

Ji(L)-i5;2(^) = e-i^^^e^(^^^-^^^+^^^?^-^^^)TV3 /" d/te-i'^^^^^^e"^^'^ 

X exp[0 (T^fc^T^/^ + ^s.PT'/^ + ^^PT'/' + k^s.T^'^ + fc^T)]. (5.77) 

Since L can be made arbitrarily small, for si E B the exponent of the term 
in the integral can be written as 

-l^T'/'il + xi) - ^2k'T{l + X2), (5.78) 

where the Xi can be made as small as desired by choosing e small enough. 
After the change of variable kT^^^ = z the integral becomes 

d^e-i^^(i+^^)e-i^^'('+^^\ (5.79) 

The integration is taken along a contour, symmetric with respect to the 
imaginary axis and such that for Re(z) > 0, arg(2) G [—n/G, 0]. This implies 
that the integral is uniformly bounded. 

Replacing the term in front of the integral ()5.77|) by one, the error can be 
estimated as 

(^e^(7^T-/3+,3,.^-V3) _ 

since the integral in ()5.77|1 is bounded. For L < e, 



(ri 
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As a consequence 

< O (r-i/Vf ('■-^)''') . (5.82) 

After this step we can also remove the error inside the integral ()5.77|) . As in 
the case of L G i3, the removal of this error leads to an additional error of 
T~^/^ with the prefactor e~i^^^^~^^^''^. Consequently we have obtained 



+ O (T-i/3e-^^('--^)''') + O (e-^^e-i^('--^)''') . (5.83) 



Next we change to the variable z = w + i(3\^ri — L. The integral becomes 



e W ' ' dw. (5.84) 

Finally completing the contour of the integration such that it goes to 
infinity in the directions arg(w) = (f± with y9+ = — vr/G and <y9_ = — 571/6 
leads to an exponentially small error. Using that 2^3 = 1//?^, the main term 
goes to ^ Ai { ^^~^ )- Since the errors are integrable in L and go to zero as 
T ^ oo, we obtain, for S2 > si, 

hm Kr{r„s,-n,s^) = f ei^(---)Ai ( Ai ^ 

T^oo \ K J \ K J 

(5.85) 

with K = \/2a3. □ 
With the change of variable X = L/n, ()5.85|1 is rewritten as 

/o 

= (5.86) 

□ 



6 The border step, Airy process 



As explained in p!lj, the Airy field is a random field which is concentrated 
on line ensembles {he{t),t G M, £ G Z+j with the properties 
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i) t hf {t) is continuous, 

ii) h^{t) < hi_i{t) for all t. 

The first line, hQ{t), of the Airy field is by definition the Airy process, denoted 
by A{t). A{t) is almost surely continuous, stationary in t, and invariant 
under time-reversal. Its single time distribution is given by the Tracy- Widom 
distribution ^Hj, known from the largest eigenvalue of GUE random matrices. 
In particular, for fixed t, 

P iA{t) >y) - e-y'^'''/^ for y ^ oo, 

F{A{t)<y) ~ e-l^l'/^2 for ?/ -oo. (6.1) 

The Airy process is localized. On the other hand it has long range correlations 
as 

{A{0)A{t)) - (A(0))2 - for \t\ oo. (6.2) 

The convergence of r]^^'^ to the Airy field, as stated in Theorem 15. H 
implies that the border step statistics, properly scaled as At, cf. p.l3j) . 
converges to the Airy process. 

Theorem 6.1. Let At{s) be the border step rescaled as in il.l!^} and let 

A{s) be the Airy process. Then for any m, Sj, G M, i = 1, . . . , m, the limit 

^lim Pt I []{AT{si) < a,} ) = P I f|{A(sA/2) < a./zt} ) (6.3) 



\i=l / V=l 



holds. 

Proof: Let /j be the indicator function of (aj, oo). Then ()6.3|) corresponds 
to 

^lim Pt {f]{VT'%f^,s^) = 0} ] = P [{\{r^^-"y{fjK,s.K/2) = 0} ) . (6.4) 



\i=l / \i=l 



We choose a large enough and split fi = f^+g"' with the indicator function 
of (oj, a] and g"- the one of (a, oo). Then 



\i=l / V=i 

m 

< 

1=1 



^T{^\{VT'VuS^)=0}) - FT[r\{VT'V-,S,)=0} 

\i=l / 

m 

Y,PT{vT"'i9'',s,)>l). (6.5) 
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The term 

PT(^nW^(/r,^0 = 0}j (6.6) 

converges to 

P (^n{^^^'"^(/r/'^, W2) = 0} j (6.7) 

which yields the right hand side of ()6.3p as a ^ oo. 

The terms in the sum of the right hand side of ()6.5p are bounded by 

Pt [vT'^a", s,) >i)<Et ^.)) = f {vr'^r, s,)) dr. 

(6.8) 

From dEIHl), 

Et {f^'\r, .,)) ^ (6-9) 

Ji{—L) is indeed a function of r + L, which asymptotics has been studied 
already for r + L large, but bounded by r + L < eT"^^'^, with the result ()5.8Hj) . 
Therefore the integrals in ()6.8|) . ()6.9|) converge for r + L < £T2/3. 

Next consider r + L > eT^/^. Let L = {r + L)T"^/^. With the change of 
variable u = ip{k) and integrating twice by parts, we obtain 



(6.10) 



Similarly as for ()5.69jl we have. 



max 

fee[i/'{-7r),i/){7r)] 



<GiL-2, (6.11) 



du'^ 

for a suitable constant Gi < oo, which yields 

\Jii-L)\<G2ir + L)-^T-'/^ (6.12) 

for some constant G2 < 00. Therefore the integrals in ()6.8j) . ()6.9|1 have a 
bound G{a) uniform in T which vanishes as a — > cxd 

□ 

Probabilistically, it would be natural to lift Theorem 16.11 to the weak 
convergence of path measures. The missing piece is the tightness for the 
sequence of stochastic process At{s). We have not attempted to fill this gap. 
The interested reader is referred to 0, where tightness for the edge scaling 
of the Aztec diamond is proved. 
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A Fermionic correlations 



A.l Two-point function 

Let A = Ak^ialai be the second quantization of the one-particle matrix 

A. It is assumed that is trace class and Det (l + e^) ^ (see [T3], Chap. 
XllI). We use the identities 



E 



E 



(A.l) 



Then 



{a*aj) 



and 



i Tr (e-\*a,) = ^ i Tr 



Z V ^ V 



e a, 



(A.2) 



nez 



Y,{a:[l + e-X,a,) = [e-%,. (A.3) 

Finally multiplying this expression by J2iez + ^~'^)~^]i,rn we obtain 

(a>,-) = [(l + e^)-i],-™. (A.4) 



A.2 Proof of (E33I)-(EI1S]) 

We prove recursively that 

■ ■ • <«in) = Det {R{ik,ji))i<k,i<n 

where 



(A.5) 
(A.6) 



Then, taking ik = jk for all k, the result ()2.43|1 - ()2.45|1 is obtained. For 
n = 1 the formula holds by definition. Suppose the formula ()A.5jl has been 
established for some n, i.e. 



(a,-, a 



(A.7) 
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We will need one more expression for (■ ■ ■ ) such that in the first k pairs the 
annihilation operator precedes the creation operator, 



/n ■ n* . . 




n* n ■ \ — 

In Jn 1 
















= (-1)'^ 
















— («i a,* ) 





(A.8) 

Let us prove this formula. For = 0, it agrees with ()A.7|) . Suppose it to be 
true for some k. Let us then prove that the formula ()A.8|) holds for k + 1, 

X^jl'^il ' ' ' '^jk + l^ik+l^ik + 2^h + 2 ■ ■ ■ '^in'^jn/ " 

= -(«J1< ■ ■ ■ + l • • • <«in) (A.9) 

+ "ik + lJk + l\^jl^il ' ' ' ^jk^ik^ik + 2^^k + 2 ■ ■ ■ ^in^jn/- 

Using the expression ()A.8j) and considering the expansion of the determinant 
in the {k + 1)*'^ column (or row), it is easy to see that ()A.9jl corresponds, 
up to a factor of —1, to the expression ()A.8|1 but with the diagonal term 
'^ife+i'^ife+i replaced by —ajk+i^^i^+i- Therefore ()A.8|) holds for k + 1, too. 

Now we prove ()A.7|) for n + 1 by using ()A.7|) for n and ()A.8|) for n and 
k < n, 

{a*a.„ ■ ■ ■ a* a,- , , ) = — Tr fe~"^a*a,-, ■ ■ ■ a* a,- , ,^ 

~ ]m,<jTr ^ Cij-^ ■ ■ ■ Oj^^^^Oj^+iO-m) 

= -Yl ("^^ii ■ ■ ■ <+i«in+i) (A.IO) 

meZ 
n+1 

+ > fe^^^l,- „ (a-., a* ■ ■ - a,- ,a* a* a,,- , , . . . a* a,-. 

p=2 

We take the term with the sum over m G Z together with the first one and 
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multiply the whole expression by J2qez + e ^jq^i^ to obtain 

(a* a,-, ■ ■ ■ a* a,- , , ) = (a* a,, ) (a* a^^ ■ ■ ■ a* a, , , ) (A. 11) 

n+1 

+ > (a* tti ) (a.- a* ■ ■ ■ a,- , a* a* a, , , . . . a* a,- , , ). 

p=2 

Using ()A.7jl and ()A.8jl for n terms we see that this last expression is noth- 
ing else than the expansion with respect to the first row of ()A.7|) with n 
substituted by n + 1. 
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